Noise-induced breakdown of coherent collective motion in swarms 
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We consider swarms formed by populations of self-propelled particles with attractive long-range 
interactions. These swarms represent multistable dynamical systems and can be found either in 
coherent traveling states or in an incoherent oscillatory state where translational motion of the 
entire swarm is absent. Under increasing the noise intensity, the coherent traveling state of the 
swarms is destroyed and an abrupt transition to the oscillatory state takes place. 
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There is a large class of problems where individual 
interacting particles, that constitute a system, are ca- 
pable of active motion and form collectively traveling 
populations. Self-propulsion of particles is already pos- 
sible in simple physical systems (see e.g. |j]-|g]) and is 
widely found in biology where individual animals may 
group themselves into swarms, fish schools, bird flocks or 
traveling cell populations |6|-[|. The role of individual 
self-propelled "particles" can also be played by localized 
patterns (spots) in reaction-diffusion systems. A bifurca- 
tion leading to the onset of translational motion of spots 
has been studied in an activator-inhibitor system with 



global feedback |l. 
diffusion systems Ml 



and in three-component reaction- 
_2[ . Interactions between individual 
self-propelled spots have been determined from the un- 
derlying reaction-diffusion equations and used to describe 
formation of bound states of such "particles" [Of . 

Mathematical modeling of collective active motion fol- 
lows several different directions. One approach is based 
on the notion of discrete stochastic automata |p|,pd]dl6[ . 
Another approach is formulated in terms of continuous 
velocity and density fields and essentially treates a swarm 
as an active fluid |L7| (such hydrodynamical equations 
may be derived by averaging from the respective au- 
tomata models fl8|]). A similar hydrodynamic approach 
is also used in the theory of traffic flows MM]]. Al- 
ternatively, one can specify dynamical equations of mo- 
tion for all individual particles that explicitely include in- 
teractions between them and/or action of external fields 
H 7||Jl^,^l| . An interesting problem related to statisti- 
cal mechanics of large populations of self-propelled parti- 
cles is the spontaneous development of coherent collective 
motion in such systems. This problem has recently been 
discussed in the framework of continuous hydrodynami- 
cal and discrete automata models, and the properties of 
the respective kinetic phase t ransiti on were numerically 

Both in one- 



and analytically investigated |17 



U 



and 



two-dimensional systems, first- and second-order transi- 
tions have been found [la . 

In the present paper we consider a population of iden- 
tical self-propelled particles near a transition between 



disordered oscillating motion and coherent translational 
motion. The particles interact via an isotropic attrac- 
tive binary potential and are subject to the action of 
noises. This globally coupled population forms a cloud 
(the swarm) in the considered one-dimensional space. 
The swarm can be found in different states. Coherent 
compact traveling states are characterized by a narrow 
distribution of velocities around a certain mean drift ve- 
locity, directed either to the left or to the right. Another 
possible state of this population corresponds to the ab- 
sence of coherent translational motion, with noisy oscil- 
lations around a certain mean position in space, deter- 
mined by the initial conditions. 

The coherent traveling states exist only for sufficiently 
weak noise and, as the noise intensity increases, the 
swarm undergoes a transition to the incoherent oscilla- 
tory state. We find that the breakdown of coherent col- 
lective motion in this system is abrupt and characterized 
by a strong hysteresis. Thus, the globally coupled swarm 
represents a multistable system that may be found in 
different states depending on the initial conditions. This 
behavior, revealed by numerical simulations, is well re- 
produced by an approximate analytical theory and may 
represent a typical property of swarms with long-range 
interactions. 

To formulate the model, we note that if a system is 
close to the onset of active motion and this instability is 
soft, i.e. characterized by a supercritical bifurcation, the 
motion with small velocity V can generally be described 
by equation 



V = aV-/3V 3 , 



(1) 



with real coefficients a and [3 > 0. This equation may be 
viewed as a normal form of the supercritical bifurcation 
leading to translational motion. Such bifurcations are 
possible in simple physico-chemical systems 0. They 
are also known for localized spot patterns in reaction- 
diffusion models and corre spond to the onset of their 
translational motion 

According to Eq. (I 
the bifurcation point (i.e 




the velocity V is zero below 
for a < 0). Above this point, 



active motion with V = ±y/a/@ is asymptotically es- 
tablished. The direction of this motion for an individual 
particle remains arbitrary and is determined by initial 
conditions. Rescaling time and introducing the new ve- 
locity variable u = V \J ' (3/a, Eq. (fil) can be written as 

u = u— u 3 (2) 

When a population of identical self-moving particles is 
considered, the velocity Ui = Xi of each particle i will 
satisfy this dynamical equation. 

Interactions between individuals may generally depend 
on both their relative positions and velocities. In this 
paper we assume that the interactions are pairwise and 
described by forces f(xi—Xj) that depend only on the dif- 
ference of coordinates of two particles i and j. We shall 
further assume that the interactions are attractive and 
depend linearly on the distance between the particles, 
i.e. f(xi — Xj) oc [xi — Xj). These attractive forces are 
supposed to model the interaction within the size ranges 
of the dynamical states considered below, where the pop- 
ulation forms clouds of either oscillating or translational 
motion. The interaction could be extended to larger dis- 
tances in order to represent, for instance, vanishing forces 
at infinity p2y . Additionally, the system may include 
noise that will be modelled by independent random forces 
£i(£) acting on individual particles. Noise prevents the 
collapse of the population, so that short range repulsion 
||15| , p2| can here be ignored. 

Under these conditions, the dynamical equations for a 
set of N identical self-moving particles with coordinates 
Xi(t) are 



JV 

Xi + {Xf ~ l)Xi + — ^ (Xi - Xj) = £i(t), 
3=1 



(3) 



for i = 1, ...JV. The coefficient a characterizes the 
intensity of interactions and can be viewed as the pa- 
rameter, specifying the strength of coupling in the pop- 
ulation. Equations (Q) constitute the basic model in- 
vestigated in this paper. We shall assume that &(£) 
are independent white noises of intensity S, so that 
(&(i)€j(t)) = 2S5ijS(t-f). Note that Eqs. (§) are 
invariant with respect to an arbitrary translation in the 
coordinate space. 

The model (||) can behave as a system of globally cou- 
pled limit-cycle oscillators (cf. p3LS3). Introducing the 
average coordinate x(t) of the swarm, 



x(t) 



1 N 
]yl>i(*). 



3=1 



Eqs. (g|) in absence of noise read 

Xi + (if — l)x l + a (xi — x) = 



(t = l,. 



(4) 



,JV). (5) 



Thus, if the swarm does not move as a whole, i.e. 
x{t) — constant, the particles perform persistent oscil- 
lations. In this state the phases of individual oscillations 
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FIG. 1. Three snapshots of 100-particle systems with 
a = 10 and S — 0.1, in different dynamical regimes. The 
central ensemble corresponds to disordered oscillations along 
a noisy limit cycle. The other two ensembles stand for coher- 
ent clouds with opposite average velocities. 



are random. Note that the spatial location x of an oscil- 
lating swarm is arbitrary. 

In addition to the random oscillatory state, the system 
(H) has two coherent collapsed states where the coordi- 
nates of all particles are identical, i.e. Xi = x for any 
i. These states correspond to uniform translational mo- 
tion of the entire swarm with the velocity u = ±1. A 
simple analysis shows that the oscillatory state and both 
coherent traveling states are linearly stable for any pos- 
itive parameter a. The final state of the population is 
determined by the initial conditions. Our numerical sim- 
ulations show that, if the average velocity u = N^ 1 J^i u i 
is initially close to zero, the oscillatory standing state is 
asymptotically reached. If, however, this initial average 
velocity is large enough, one of the two coherent traveling 
states will be approached. 

Since the particles either converge to coherent motion 
with constant velocity or to disordered oscillations with 
no average drift, the ensemble can be thought of as a 
multistable system with qualitatively different attractors. 
In the following, we focus the attention on how these at- 
tractors respond to the effect of noise. With this aim, 
we study Eq. (ph numerically. Integration is performed 
by means of a standard Euler scheme with a time step 
At = 1CP 3 to 10~ 2 . Most calculations correspond to en- 
sembles of 100 particles, with the coupling intensity rang- 
ing from a = 1 to 100. Larger values of a require smaller 
values of At. Noise is introduced by generating, at each 
time step, a random number £ with uniform distribution 
in the interval (— £o> £o)- This choice corresponds to hav- 
ing S = £,o/6At. In practice, £ is calculated for each 
given value of S. Initial conditions are selected at ran- 
dom, distributing the particles around x = and u = 
or 1 with a dispersion of the order of 0.5 in both vari- 
ables. From each initial condition the system is left to 
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FIG. 2. Normalized profiles of coherent clouds as functions 
of the deviation from the average position, for different values 
of a and S — 0.1 in a 100-particle ensemble. 



evolve in the absence of noise until it reaches the state 
of disordered oscillations or coherent motion. Then, at 
t = 30, noise is switched on. Typical calculations extend 
up to t w 1000. 

For small noise intensities, S < 0.1, the stochastic per- 
turbations to the trajectories preserve the characteristic 
features of the collective dynamics observed in the ab- 
sence of noise. The completely collapsed state of the 
noiseless case transforms into a cloud of particles which 
still moves coherently at a given velocity. Oscillatory 
orbits, meanwhile, proceed now along a noisy limit cy- 
cle. Figure 1 shows three snapshots of a system of 100 
particles with a = 10, subject to noise with S = 0.1. 
They started from different initial conditions, as de- 
scribed above. The arrows indicate the overall motion 
of each swarm. 

Within coherent clouds, each particle performs an os- 
cillatory noisy motion which is superimposed to the col- 
lective translation. The distribution of particles inside 
the clouds has a well defined profile, shown in Fig. 2 
for some values of a in the case of positive velocity. The 
normalized distribution p(y) is there plotted as a func- 
tion of the coordinate relative to the average position, 
Hi = Xi — x. For decreasing a, the distribution becomes 
broader and more asymmetric, with an accumulation of 
particles at the front of the cloud. 

The coherent traveling states of the population cease to 
exist at sufficiently high noise intensities and the swarm 
undergoes an abrupt transition to its random oscillatory 
state, characterized by the absence of the translational 
motion. This breakdown of coherent swarm motion is 
illustrated in Fig. 3. We see that if the noise is rela- 
tively weak (Fig. 3a), switching it on at t = 30 only 
produces a slight decrease of the velocity of the coherent 
cloud, so that the average velocity u(t) exhibits fluctua- 
tions around a constant mean value u < 1. If however 



FIG. 3. Average velocity of 100-particle coherent clouds in 
two realizations at (a) S = 0.10 and (b) 0.12, with a = 10. 
Noise is switched on at t — 30 (dashed line). The horizontal 
lines indicate the asymptotic mean values of u(i). 



the noise intensity exceeds a certain threshold, the effect 
of introducing noise is qualitatively different (Fig. 3b). 
Within a certain time interval after the introduction of 
noise, the swarm continues to travel at a somewhat re- 
duced, strongly fluctuating average velocity u(t). Then, 
it suddenly starts to decelerate and soon reaches a steady 
state where the average velocity u(t) fluctuates near zero. 
Inspection of the distribution of particles in the ensemble 
shows that in this state the system has been attracted to 
the noisy limit cycle mentioned above. We conclude that 
the system undergoes a noise-induced transition from a 
condition of multistability with two kinds of attractors to 
a situation where only one of them exists. The coherent 
clouds observed for small noise intensities are no longer 
possible for S > S Cl and the system is necessarily led to 
the state of noisy, disordered oscillations. 

Figure 4 displays the dependence of the mean velocity 
u of the traveling swarm on the noise intensity S for three 
different values of the coupling coefficient a. We see that 
the mean velocity monotonously decreases with the noise 
intensity, until a certain critical noise intensity is reached 
and the coherent swarm motion becomes impossible. The 
mean velocity at the critical point is still relatively large, 
u ~ 0.8. The critical noise intensity S c becomes lower for 
smaller values of a. Note that the behavior of the swarm 
is characterized by a strong hysteresis. If the breakdown 
of the coherent motion has occurred, subsequently de- 
creasing the noise intensity leaves the system in the os- 
cillatory state with zero mean velocity, down to S = 0. 

An interesting property of the considered noise- 
induced transition is the divergence of the waiting time 
at the critical point. The waiting time To is defined as 
the time at which the average velocity u(t) of the cloud 
first reaches zero (we measure this time starting from the 
moment t = 30 when the noise is switched on). Figure 5 
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FIG. 4. The asymptotic mean velocity of 100-particle co- 
herent clouds as a function of S, for different values of a. Dots 
correspond to numerical measurements and lines stand for the 
analytical result Eq. did). 



FIG. 5. Waiting time To, at which the average velocity 
of initially coherent clouds vanishes for the first time, as a 
function of S — S c for a — 10 in a 100-particle system. The 
dashed lines have slopes —1/2 and —4/3. 



shows the waiting time To as a function of S— S c in a log- 
log plot. We see that for very small values of S — S c , this 
time decreases following a power law, T ex (S — S^) -7 , 
with 7 rs 1.33. Then, at about S — S c = 0.03, the be- 
havior changes to a power law with 7 rs 0.52. Straight 
dashed lines with slopes —4/3 and —1/2 have been plot- 
ted for reference. 

The observed noise-induced transition between coher- 
ent clouds and disordered oscillations of the swarm can be 
explained by a simple approximate analytical approach. 
By summing all Eqs. (o) for different particles i and 
taking into account that the noises acting on individual 
particles are not correlated, an evolution equation for the 
average swarm velocity u{t) is obtained: 



1 



N 



N *—> 



.3 T7 



U = 0. 



(0) 



Let us introduce for each particle its deviation y+ = Xi — x 
from the average position of the swarm. Then we can 
write 



1 N 



N 



3era 



1 



i=i 



N 



(7) 



where a = TV -1 J2i vl i s the average square dispersion of 
the swarm. The last cubic term in this equation can be 
neglected if the distribution of particles in the traveling 
cloud is symmetric. As we have seen from numerical 
simulations (Fig. 2), this is indeed a good approximation 
for sufficiently large values of the coupling constant a. 
Within this approximation, Eq. (0) takes the form 



u + (u — l)u + 3au = 0. 



(8) 



On the other hand, deviations of particles from the cen- 
ter of the swarm obey the stochastic differential equation 



jji + (3u 2 - l)jji + ay t + 3m (y. 




&(*). 



(9) 



Assuming that the deviations of yi are relatively small 
and linearizing this equation, we obtain 



jji + (3m 2 -l)yi + ayi = &(i). 



(10) 



In this approximation the deviations for different par- 
ticles i represent statistically independent random pro- 
cesses. This allows us to replace the ensemble average in 
the dispersion a by the statistical average taken over in- 
dependent random realizations of such processes, defined 
by Eq. ©. 

Hence, we have derived a closed set of equations <M) 
and (|10|) that approximately describe the swarm. We 
want to investigate steady statistical states of this sys- 
tem. The stationary solutions to Eq. (||) are u = 
±\/l — 3cr and u = 0. The latter solution corresponds 
to the resting swarm. 

Examining Eq. (QOh , we note that it describes damped 
oscillations only if 3lP — 1 > 0, i.e. only if mean velocity 
of the swarm is sufficiently large. Under this condition, 
the stationary probability distribution for y^ is readily 
found and the average square dispersion of velocities is 
obtained as 



(11) 



3m^ 



1 



The algebraic equations for u and a can be solved, 
yielding the statistical dispersion of particles in the trav- 
eling swarm, 



a l>2 = i (l ± %/T - IS) 



(12) 



and its mean velocity 



—2 



n,2 



\(2±VT 



9S 



(13) 



Thus, the traveling state solutions disappear when the 
critical noise intensity S c = 1/9 = 0.11 ... is reached. 
At this critical point the mean swarm velocity is u c — 
a/2/3 = 0.82 . . . and the mean dispersion of particles in 
the cloud is <j c = 1/9 = 0.11 . . .. 

Below the breakdown threshold (for S < S c ), solution 
(O) has two branches shown by solid and dashed lines in 
Fig. 4. The lower branch is apparently unstable, since it 
approaches the value u = l/-\/3 = 0.58 . . . at S = 0, i.e. 
in absence of the noise. A special property of the derived 
solution is that it does not depend on the parameter a. 

Comparing the theoretical prediction with the numeri- 
cally determined values of the mean swarm velocity, that 
are also plotted in Fig. 4, we can see that this approx- 
imation provides good estimates of the swarm velocity 
and the critical noise intensity when the parameter a is 
relatively high (a = 100 and a — 10). At small values 
of a, the deviations from the numerical results become 
significant near the breakdown threshold. This can be 
understood if we take into account that, according to 
Fig. 2, the distribution of particles in a traveling swarm 
shows significant asymmetry for such a small value of a 
and therefore our approximations are not valid. 

For a standing swarm (u = 0), the deviations yi = 
Xi — x obey in the limit N —* oo the nonlinear stochastic 
differential equation 



Vi + (yf - 1)& + a Vi = &(*), 



(14) 



which is similar to the Van der Pol equation |25j . In this 
state, therefore, the particles in the swarm perform pe- 
riodic limit-cycle oscillations with a random distribution 
of phases. This state exists for any noise intensity S and 
is approached when the noise-induced breakdown of the 
coherent motion takes place at S = S c - 

Thus, we have find in this paper that a swarm of inter- 
acting, actively moving particles may show bistable be- 
haviour, i.e. be found either in a coherent state, traveling 
at a fixed velocity, or in a rest state where the transla- 
tional motion is absent and the individual particles per- 
form oscillations around the center of the swarm. The 
bistability persists in the presence of noise if its inten- 
sity remains relatively low. Increasing the noise intensity 
leads to a sudden breakdown of the coherent traveling 
motion and a transition to the resting oscillatory state 
occurs. This behavior is different from the second-order 
phase transitions to coherent collective motion, that were 
found in the previously studied models |17|Jlq| . We con- 
jecture that the difference is related to the fact that in 
our model the interactions between self-propelled parti- 
cles have a long range and extend over the entire swarm. 
It would be interesting to see how this behavior is modi- 
fied when other interaction laws and systems with higher 
dimensionality are considered. Finally, we remark that, 



when formulated in terms of dynamical equations for in- 
dividual interacting self-propelled particles, the problem 
shows significant similarities to synchronization and con- 
densation in populations of globally coupled oscillators 
(see e.g. J23LE6|). The significant new aspect is that 
collapsed synchronous states correspond here to trans- 
lational motion of the entire population. 
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